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Abstract 

Context. Imaging polarimetry is a powerful tool for imaging faint circumstellar material. It is a rapidly developing field with great 
promise for diagnostics of both the large-scale structures and the small-scale details of the scattering particles. 
Aims. For a correct analysis of observations we need to fully understand the eff'ects of dust particle parameters, as well as the eff'ects 
of the telescope, atmospheric seeing, and assumptions about the data reduction and processing of the observed signal. Here we study 
the major eff'ects of dust particle structure, size-dependent grain settling, and instrumental properties. 

Methods. We performed radiative transfer modeling using difi'erent dust particle models and disk structures. To study the influence 
of seeing and telescope diff'raction we ran the models through an instrument simulator for the ExPo dual-beam imaging polarimeter 
mounted at the 4.2m William Herschel Telescope (WHT). 

Results. Particle shape and size have a strong influence on the brightness and detectability of the disks. In the simulated observations, 
the central resolution element also contains contributions from the inner regions of the protoplanetary disk besides the unpolarized 
central star. This causes the central resolution element to be polarized, making simple corrections for instrumental polarization diffi- 
cult. This eff'ect strongly depends on the spatial resolution, so adaptive optics systems are needed for proper polarization calibration. 
Conclusions. We find that the commonly employed homogeneous sphere model gives results that differ significantly from more re- 
alistic models. For a proper analysis of the wealth of data available now or in the near future, one must properly take the eff'ects of 
particle types and disk structure into account. The observed signal depends strongly on the properties of these more realistic models, 
thus providing a potentially powerful diagnostic. We conclude that it is important to correctly understand telescope depolarization and 
calibration effects for a correct interpretation of the degree of polarization. 
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1. Introduction 

High spatial-resolution imaging polarimetry is a rapidly devel- 
oping field. It is now recognized that imaging polarimetry is 
a powerful tool for imaging faint matter around a bright star. 
Polarimetric diflferential imaging (PDI) is used to suppress the 
glare of the central star, which is unpolarized, to be able to see 
the faint reflected light from the material surrounding the star, 
which becomes polarized by scattering. This has already been 
successfully used to image circumstellar disks around young 
stars (see e.g. [Close et al.| 1997] [Hashimoto et aL 



2011) and cir- 



cumstellar dust shells around evolved stars (e.g [Gledhill et al. 
|2Q01[ ), and it has been proposed for imaging the light reflected 
from exoplanets ( Stam et al.|2Q04| ). Images in polarized intensity 
not only contain information on the geometry of the circumstel- 
lar material, but also carry information on the particles scatter- 
ing the light from the central star. In the general case, these two 
characteristics are intertwined and not easily separated. 

Perhaps the most challenging case for interpreting polarimet- 
ric images is that of optically thick circumstellar disks around 
newly born stars. First, being optically thick causes multiple 
scattering, which alters the degree of polarization. The geom- 
etry of the surface of these disks can also cause scattering at 
unex pected angles , in turn causing certain features in polariza- 
tion ( Perrin et al.[[20b9| ). These eflTects can complicate the in- 
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terpretation of polarimetric images in terms of planet-disk in- 
teractions. Besides these geometrical eflTects on a large scale, a 
second complicating factor arises when these images are used to 
derive general information about the dust particles in the disk. 
Dust particles are expected to be highly complex in shape and 
might be quite large in size, since they are a mixture of smal l 
grains and large, complex fluffy aggregates ('Dominik"2009[). 
Computing the optical properties of such particles is very dif- 
ficult and computationally demanding, if feasible at all. In ad- 
dition, the diflTerent grain sizes are not expected to be mixed 
well with the gas, but size-dependent settling of the dust par- 
ticles causes a size distribution of the grains that changes with 
height in the disk ( [Dullemond & Dominik||20Q4[ ). To conclude 
this chain of complications, dust properties, i.e. gr ain size and 
composition, are expected to vary radially (see e.g. van Boekel 
et al. 2004). Previous studies have already covered large parts 
of the parameter space of dusty circumstellar disks in the con- 
text of polarimetric imagin g. The eflTects of disk mass and radial 
extend were studied in e.g. [Whitney & Hartmann[ ( p^992| ). More 
recently, Murakawa (2010 ) has shown computations of near in- 
frared polarimetric imaging of disks with various grain sizes and 
parameterized thickness of the disk. 

An additional complication is the eflTect of the finite reso- 
lution of real telescopes on the final image. The most power- 
ful diagnostic, the degree of polarization of the scattered light, 
is highly influenced by the glare from the central star. Also, 
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blurring of the high-intensity, polarized, inner regions of the 
disk will affect the appearance of the disk in polarized inten- 
sity. Furthermore, data reduction techniques that correct for in- 
strumental polarization using the central star as an unpolarized 
calibration point might suffer from the non-zero innermost disk 
polarization. 

In this paper we present computations of the polarimetric im- 
ages of protoplanetary disks around young, intermediate-mass 
stars, so-called Herbig Ae stars. We study how the images are 
effected by the structure of the dust grains and size-dependent 
grain settling. We varied the size of the dust particles, the geom- 
etry of the disk, and the shape of the dust particles. This was done 
using state-of-the-art models for the optical properties of the dust 
particles combined with a highly efficient radiative transfer code. 
An instrument simulator for a dual-beam imaging polarimeter at 
a 4 m class telescope was used to study the effects of atmospheric 
seeing and data reduction artifacts. 

2. Model setup and computational approach 

In this section we describe the setup of the models in terms of 
dust grains and disk geometry, in Section |2.1| we motivate our 



choice for the parameters used in the computations. The radia- 
tive transfer model used to create the polarimetric images and 
the instrument simulator are also described. 



2.1. Parameter space 

The characteristics of dusty, protoplanetary disks are described 
by a large set of parameters ranging from the microscopic prop- 
erties of the dust grains to the macroscopic g eometrical prop- 
erties of the disk. [Whitney & Hartmann ( 1992 ) investigated the 
effects of disk mass, radial extent, and the amount of disk flar- 
ing. They conclude that a flat disk would not be detectable. Also, 
they conclude that the large extent of many observations avail- 
able at that time was most likely caused by an extended infalling 
envelope. The appearance of such an envelo pe in combination 
with the effects of outflows was studied in Whitney & Hartmannl 
( |1993| ) and [Fischer et al.| ( ^1996 ) for a wide range of enve lope 
and outflow parameters. More recently, Murakawa ( 2010| ) has 
presented computations of near infrared polarimetric imaging of 
disks. These theoretical studies all provide a good analysis of the 
dependence of polarimetric imaging on certain parameters of the 
dusty surroundings of young stars. It is therefore not the aim of 
this paper to repeat this and present a study covering the entire 
parameter space. 

There are a few aspects that have not received much attention 
in previous studies, so we focus on these aspects and study their 
effects on the observed polarimetric images. These aspects are 



- Grain structure, 

- Size dependent grain settling, 

- Atmospheric blurring, and other observational aspects. 

Below we explain how we set the model up specifically keeping 
these three main points in mind. 

2.2. Dust grain model 
2.2.1 . Composition 

The composition of protoplanetary dust is far from trivial (for 
a recent review see Min & Fly nn 2010 ). To arrive at a simple, 
yet quite accurate model of the com position, we follow the pro- 
cedure also adopted by Min et al. ( [2011i ). First we define the 



abundances of all available dust-for ming elements to be the sam e 
as the solar composition derived by'Grevesse & Sauval ( 1998 ). 
In meteorites and interstellar dust particles (IDPs) in our solar 
system all of the available sulfur is in the form of iron sulfide. 
Therefore, we first put all S in FeS thereby taking away about 
half of the available Fe atoms. From the remaining Fe, Mg, and 
Si we produce an amorphous Fe/Mg silicate with a stoichiome- 
try in between olivine and pyroxene. We model this mixture by 
averaging the optical properties of different types of amorphous 
silicates. In addition to this iron sulfide/silicate mixture we add a 
contribution from carbonaceous matter. It is found that in our so- 
lar system about half of the carbon atoms are in the solid phase 
(Geiss 1987 J , so we add this abundance in the form of amor- 
phous carbon grains. The composition resulting from this con- 
tains by mass 58% silicates, 18% iron sulfide, and 24% amor- 
phous carbon. 

To compute the opacities of the dust grains, we need the re- 
fractive index as a function of wavelength. This has been mea- 
sured in t he laboratory. For the s ilicates we use the measure- 
ments by |Dorschner et al. (11995') and |Henning & Stognienko 
(199 6|), for the iron su lfide we take the laboratory data from 
Beg emann et al.| ( [1994 ), and for the carboneceous material we 
take the measurements from Preibisch et al.| ( p^93] ). 

It can be argued that grains in the outer regions of protoplan- 
etary disks contain a coating of ice. However, there is a recent 
debate about whether ice layers can survive in the upper layers 
of the di sk where it is exposed to UV radiation from the central 
star (seel Grigorieva et al.|2007| ). We checked how the scattering 
properties are influenced when adding water ice to the mixture in 
an abundance of 45% by mass, which is the abundance expected 
when putting all remaining oxygen in water ice ( Mi n et al. 2011| ). 
The scattering properties are affected only slightly. The single 
scattering albedo goes up by less than 15% in all cases. The ex- 
tinction coefficient, the maximum degree of polarization, and the 
asymmetry of the phase function are hardly affected. The result- 
ing polarimetric images are thus very similar to those without 
ice. 

2.2.2. Grain structure 



We distinguish between two extreme types of particle structures 
using the above composition: 

1. compact, homogeneous grains, and 

2. fluffy, inhomogeneous aggregates. 

These two different particle structures require different ways of 
computing the optical properties. To correctly estimate the scat- 
tered light intensity and degree of polarization, the model used 
for the optical properties of the grains is crucial. Homogeneous 
spherical grains cause resonances at certain scattering angles and 
generally behave differently from more natural, irregular shapes. 
However, for reference and comparison to previous work, we 
also computed several models using compact, homogeneous 
spheres. 

Compact, homogeneous grains: This model considers a mix- 
ture of grains, each composed of a single material. The optical 
properties of such a mixture can be obtained by simply averag- 
ing the optical properties of homogeneous grains using the given 
abundances of the various materials. We chose here to model 
the optical properties of the dust grains using the statistical ap- 
proach. In this approach we take a wide distribution of simple 
particle shapes to simulate the variety of irregularities occurring 
in nature. We use the distribution of hollow spheres (DHS; see 
Min et al. 2005 1 ) with an 'irregularity parameter' /max = 0.8. 
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Figure 1. The scattering properties of the dust particles for different values of the maximum particle size, <2max. Upper panels: the 
normalized phase function, Fn, as a function of scattering angle, 0. Lower panels: the degree of linear polarization as a function of 
scattering angle. The wavelength is A = 600 nm. Solid lines are for the fluffy aggregates, dashed lines for the compact grains. 



It turns out that this choice of this grain-shape model, together 
with our choices for the grain sizes and compositions, provides 
a scattering matrix that is quite similar to cometary dust. 

Fluffy, inhomogeneous aggregates: When small dust grains 
collide, they can form large aggregated structures. The optical 
properties of such aggregates display characteristics of both the 
constituents and the size of the aggregate as a whole (see e.g. 
Volten et al.| |2007|). We use here the approximate method by 
Min et al. 72008 ) to compute the optical properties of large, 
inhomogeneous fluffy aggregates. In this method we use effec- 
tive medium theory to mix the dust material with vacuum and 
thus obtain a highly porous, fluffy structure. The vacuum frac- 
tion can be computed from the aggregate size and fractal dimen- 
sion. Though this method was not primarily designed to simulate 
scattering at optical wavelengths, we find that many of the char- 
acteristics attributed to fluffy aggregates are simulated well. In 
general we find that the grain s computed using these methods 
with the parameters given by |Min et al.| ([2008 ) are extremely 
fluffy and as such give an extremely forward-peaked scattering 
phase function and a very high degree of polarization at a 90° 
scattering angle. 

2.2.3. Grain size 

To address the effects of grain growth on the observable polari- 
metric signal we use a parameterized grain size distribution. We 
adopt a power law size distribution with n{a) oc a~^ . Here a is the 
volume equivalent radius of the grains and the distribution runs 
from amin to Omax- We fix amin to 0.03 yum and vary a^ax from 
0.3 yum to 3 mm. The different grain sizes are treated separately 
in the radiative transfer. This implies that they can have different 
temperatures and a different vertical distribution caused by grain 
settling. 

Figure[T]displays the single scattering phase function and the 
degree of linear polarization as functions of the scattering angle, 
6, for different values of ^max at a wavelength of A = 0.6 /jm. 
The most important scattering characteristics are summarized in 
Table [1] Shown are the mass extinction coefficient, /Cext, the sin- 
gle scattering albedo, oj, the asymmetry parameter of the phase 



Particle type 


^max 


^ext 




(cos 6) 


p 

^ max 




[//m] 


[cmVg] 






[%] 


Compact grains 


0.3 


35290 


0.63 


0.53 


33 




30 


1532 


0.73 


0.76 


40 




3000 


25 


0.71 


0.79 


44 


Aggregates 


0.3 


31246 


0.51 


0.87 


95 




30 


5277 


0.57 


0.98 


89 




3000 


139 


0.53 


0.99 


89 



Table 1. Light scattering characteristics of the different particle 
types for a few values of ^max. 



function, (cos 6), and the maximum degree of polarization from 
single scattering, Pmax- 

2.3. Disk model 

For the disk model we try to avoid parameterization as much as 
possible. Therefore we compute the vertical scale height of the 
disk self-consistently. This means that the vertical density dis- 
tribution is fully determined by the temperature structure of the 
disk assuming vertical hydrostatic equilibrium (see Dulle mond] 
et al. 2007 for a review of the theoretical picture of the structures 
of protoplanetary disks). In addition, we compute the radial lo- 
cation and structure of the inner rim using ev aporation and sub- 
limation physics as described in Kama et al.| ([2009 ). This results 
in an inner edge of the disk that is curved outwards and slightly 
higher than the regions behind it. For details we refer to Kama] 
^taL| ( [20091 ). 

For the central star we take a Herbig Ae star. We use two 
different disk types, a flaring disk, where we assume gas and 
dust to be perfectly coupled, and a settled disk, where we use 
the settling prescription from |Dubrulle et al.] ( |1995| ) as described 
above. These two disk types a pproximately represe nt the Group 
I and II sources as defined by 'Meeu s et al.| ( [200T] ). All star and 
disk parameters are listed in Table [2] 

As the grains grow to larger s izes they decouple from the ga s 
and settle towards the midplane ( [Dullemond & Dominik|2004 ). 
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parameter 


symbol 


value 


mass of the star 




2.5 Mo 


effective temperature 




10000 K 


radius of the star 




2Rq 


dust mass in the disk 


Mdust 


1 • 10-3 Mo 


distance to the star 


D 


ISOparsec 


outer radius of the disk 


Rout 


2000 AU 


surface density turnover point 


Ro 


200 AU 


minimum grain size 


^min 


0.03 fim 


maximum grain size 


^max 


0.3,30,3000//m 


turbulence parameter 


a 10- 


(settled), oo (flared) 


gas-to-dust ratio 


fg/d 


120 


Table 2. Values for the model 


parameters 





This is thought to be a possible explanation for the flatter disk 
structures observed around some stars. A full treatment of set- 
tling is beyond the scope of this paper and will be implemented 
in the future in the MCMax radiative transfer code. Here we pa- 
rameterize settling using the approximate method by [Dubrulle 
et al.| (p'995 ). In this way we do simulate the dependence of set- 
tling efl&ciency on grain size, grain density, gas density, and gas 
temperature. To distinguish between effects of grain settling and 
grain optical properties we take for the grain settling parameters 
always those of the material volume equivalent sphere. In prin- 
ciple, settling is more efficient for homogeneous compact grains 
and less efficient for fluffy aggregates. This is simply because 
the larger surface area of the fluffy grains allows them to cou- 
ple much more easily to the gas. However, since we want to be 
able to separate disk structure effects as much as possible from 
effects of grain properties, we have chosen to not let this effect 
influence the settling efficiency and determine the settling prop- 
erties of all grains as if they were spherical homogeneous grains 
with the same mass. The description of grain settling we use has 
one free parameter, which is the strength of the turbulence, a. 
This parameter defines the balance between settling and vertical 
mixing. For the settled disk models we take a = 10~^, while for 
the well mixed models we set a = cx). 

The radial density distribution of the dust disk was pa rame- 
terized using a radial surface density ( Hughes et al.|2008] ) 



E(r) oc r ^ exp < - 



2-p 



(1) 



law in the inner region. We fix this to p = I, wh ich is a com- 
monly used value (see e.g. DuUemond et al.|2()06| ). Although we 
take the outer radius of the disk to be 2000 AU, the surface den- 
sity of the disk is already highly diminished after the turnover 
radius Rq = 200 AU. 

In Fig. [2] we depict the two different disk structures consid- 
ered. This figure illustrates that the scattering angles of the sur- 
face layers of the disk depend strongly on the opening angle. In 
the case of a very flat disk, the scattering angles observed at the 
front and the back of the disk are symmetric around 90°. With 
/ the inclination angle of the disk the scattering angles at front 
and back for a very flat disk are 6 = 90° - / and 6 = 90° + / 
respectively. For a strong flaring disk with a large opening angle 
j/^, the scattering angles at front and back are symmetric around 
90° - i// and become ^ = 90° - - / and 6/ = 90° - + /. 

We would like to emphasize that by using the self-consistent 
treatment described above we limit the amount of free param- 
eters significantly. Besides the parameters we do vary, we have 
two important parameters that determine the observed polariza- 
tion signal; Rq and Mdust- The impact of these two parameters 
on t he polarimetric signal o f circumstellar disks was studied in, 
e.g., [Whitney & Hartmann| ( |l992| ). 



for r < /?out- Here Rq is the turnover point from where an expo- 
nential decay of the surface density sets in and p sets the power 



2.4. Radiative transfer 

The transfer of radiation through the disk is done in full using the 
3D axisymmetrical radiative transfer code MCMax ( |Min et al.| 
i2Q09j . This code uses Monte Carlo radiative transfer to compute 
the disk temperature structure and from this the vertical density 
scale -height and the sublimation state of the dust (see also'KamaJ 
|et al . 2009 ). The scattering of radiation uses the full scattering 
matrix acting on the Stokes vectors. The grains in the disk are 
assumed to be in random orientation so symmetry arguments can 
be used to simplify the scattering matrix. 

Often images and other observables are extracted from 
Monte Carlo radiative transfer methods by simply collecting 
the escaping photons after they are traced through the disk. 
Although this is an elegantly simple way, we adopt a different 
technique to construct the images which allows us to obtain ac- 
curate results even in the very faint parts of the image. After the 
temperature structure of the disk is obtained, images are con- 
structed by integrating the formal solution to the radiative trans- 
fer. For this, the local scattering source function is also needed; 
i.e., at each volume element in the disk we need to know how 
much radiation is scattered towards the observer. This is again 
computed using a Monte Carlo approach. However, this time we 
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store the full 3D radiation field at each location in the disk by 
using the full path traveled by the Monte Carlo photons (as in- 
spi red by the p ath length method for temperature determination 
by Lucy 1999| ). In combination with the local thermal source 
function, which is obtained from the temperature, the formal so- 
lution can be easily integrated to obtain smooth images. This is 
done by integrating from each volume element in the disk the 
total emitted flux (as obtained from the local temperature) and 
the total scattered flux (as obtained from the 3D radiation field) 
diminished by the extinction through the disk along the line of 
sight towards the observer. 

The advantage of the above method over 'classical Monte 
Carlo' photon-collecting methods is that the information accu- 
mulated during the Monte Carlo photon-tracing process is used 
to the fullest to obtain the local 3D radiation field. This way 
smooth images can also be produced in regions where only a 
small fraction of the photons arrive. This is especially important 
for imaging polarimetry where one of the main purposes is to 
gain a very high dynamic range in the image. The method was 
benchmarked against other radiative transfer codes in Pinte et al. 
( |2QQ9l ). 

2.5. Instrument simulation 

Ground-based observations are aff'ected by atmospheric and in- 
strumental eff'ects such as seeing and instrumental polarization. 
To account for these eff'ects, a numerical simulation includ- 
ing readout noise, photon noise, instrumental polarization, and 
diff'erent seeing conditions was performed. These simulations 
model the dual-beam imaging polarimeter ExPo (the Extreme 
Polarimeter; Rodenhuis et al., 2008^ 2011) , currently a visitor in- 
strument at the 4.2 meters William Herschel Telescope (WHT). 
ExPo measures linear polarization at optical wavelengths, and 
combines very short exposure times with the dual-beam tech- 
nique to minimize systematic errors. In a standard observation, 
thousands of single exposure frames are combined to obtain a 
total exposure time of ~ 15 minutes. 

The point spread function (PSF) for ExPo is calculated as 



PSF(^) = \r{A • e' 



(2) 



where stands for fourier transform, A and (p refer to the ampli- 
tude and phase of the incoming wavefronts, respectively. A set 
of 100 statistically independent wavefronts is calculated for dif- 
ferent seeing conditions. We compute the phase of these wave- 
fronts according to the Kolmogorov theory of turbulence (for 
details see Appendix |A.1[ ). Broadband PSFs are obtained from 
these sets of wavefronts combined with the aperture of the tele- 
scope after which they are convolved with full-resolution simu- 
lated images. Readout noise, photon noise, instrumental polar- 
ization, and guiding errors are added to these images. This set 
of simulated-raw images is then reduced with the ExPo Pipeline, 
described in detail in Canovas et al.| ( |20 1 1). 

The characteristics of this simulation are summarized in 



Table A. 1 The simulation of the atmosphere, telescope, and in- 
strument is detailed in Appendix [A[ 

3. Resulting full-resolution images 

In this section we discuss the intrinsic polarimetric images of 
the disks and the eff'ects of the diff'erent dust properties and disk 
parameters. All full-resolution model images are generated at a 
wavelength of A = 600 nm. Though detailed studies of the color 
behavior of the polarized images can be a useful diagnostic, the 



general monochromatic behavior of the images is not very sen- 
sitive to the exact wavelength. We use these single wavelength 
images to represent the average broadband signal, so the broad- 
band eff'ects are fully taken into account for the simulation of 
the telescope and instrument, but the full-resolution model disk 
images are computed at a single wavelength. 

In Fig. [3] we present polarization maps for the diff'erent dust 
types. We show here images for our two main grain types (com- 
pact grains and flufl'y aggregates) and also for compact homoge- 
neous spherical grains. It is clear from this figure that the type 
of dust grain has a strong influence on the appearance of the 
disk. Especially in the upper part of the image, the homogeneous 
spherical dust grains display a behavior not seen for the more re- 
alistic dust models. The homogeneous spherical grains display 
here a decreasing polarized intensity towards the upper central 
region of the disk. This feature is seen in a much milder form in 
the image of the flufl'y aggregates, but generally the images for 
the other two grain shape models show smoother behavior in this 
part of the image. Also, it is clear that the flufl'y aggregates are 
fainter than the compact grain models. This is caused by the ex- 
tremely forward-scattering nature of the flufl'y aggregates, where 
only a small fraction of the light is scattered towards the ob- 
server. 

In Figs. |4]and|5] we show the polarized intensity maps for 
the compact grains and the flufl'y aggregates for difl'erent val- 
ues of the maximum grain size in the size distribution. For the 
flared disk models (Fig. |4]), it is clear that the brightness and 
spatial extent of the disk decreases when the particle size is in- 
creased. This is caused by the reduced opacity of the larger parti- 
cles, which decreases the flaring angle of the disk. For the settled 
disk models we see a slightly difl'erent trend. Because the surface 
layer of the disk is here already in the shadow, the decrease in 
the height of the disk is not so important. More important in this 
case is the increase in the scattering albedo, causing the disk to 
become brighter and appear larger at first when the particle size 
is slightly increased. When the particle size is increased further, 
the decrease in the total opacity becomes dominant and the disk 
becomes fainter again. 



Murakawa] ( |20 1 0| ) revived the discussion on the 'roundabout 
eff'ect' for edge on disks in near infrared polarimetric imaging. 
In the dark lane of an edge on disk no first-order scattering can 
be observed since no starlight can reach these regions in the disk. 
Second-order scattering causes the polarization vectors not to be 
aligned perpendicular to the direction towards the central star, 
bu t parallel to the dark lane o f the disk. Th is effect was first noted 
by|Bastien & Menard] ( |1 98 8| ), after which [Whitney & Hartmann 
(1992) argued that detecting this feature would be very hard be- 
cause of the low surface brightness in this part of the image. We 
also observe this eff'ect in our optical polarization images of the 
most flaring edge-on disks. Figure [6] shows the polarized inten- 
sity maps for the disks seen edge-on for difl'erent grain sizes. We 
only show the compact grain models, but the flufl'y aggregate 
models show similar behavior. Clearly we see the roundabout 
efl'ect in the flaring disk with only small grains. As soon as the 
grain size is increased, or the disk is settled, the efl'ect disappears 
because the dark lane becomes too narrow. The region where 
the efl'ect can be observed is very weak in polarized intensity as 
compared to the region above and below. This implies that, in 
observations with a limited spatial resolution, the effect will be 
washed away by PSF smearing. 
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Figure 3. Comparison of the polarized intensity maps for different grain types. Upper panels are for the flared disks, lower panels 
for the settled disks. Left panel: homogeneous spheres; middle panels: compact, irregularly shaped grains; right panels: fluffy 
aggregates. For all images, the maximum size of the particles is 30yum. Spatial scale is in arcseconds. The length of the vectors is 
proportional to the degree of polarization. The inclination of the disk is 55° with respect to pole on. 
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Figure 4. Comparison of the polarized intensity maps for different grain sizes. All panels are for the well mixed, i.e. flaring, disk 
models. Upper panels are for the compact grains, lower panels for the fluffy aggregates. The maximum particle size in the size 
distribution increases from left to right. Spatial scale is in arcseconds. The length of the vectors is proportional to the degree of 
polarization. The inclination of the disk is 55° with respect to pole on. Color scale is the same as in Fig.|3] 
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Figure 5. Same as Fig.|4j but for the settled disk models. 
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Figure 6. Comparison of the polarized intensity maps for different grain sizes in the nearly edge on case. All panels are for the 
flared disk models with compact dust grains. The settled disks display only a very thin signature. The maximum particle size in the 
size distribution increases from left to right. Spatial scale is in arcseconds. The length of the vectors is proportional to the degree of 
polarization. The inclination of the disk is 85° with respect to pole on. Color scale is the same as in Fig.|3] 



4. Simulated observational images 

We present images simulated for the ExPo instrument mounted 
at the 4.2 m William Herschel Telescope (WHT). The aim here 
is to get a feeling for the effects of atmospheric seeing, telescope 
diffraction, and noise on the interpretation of polarimetric im- 
ages. 

The final image observed through a telescope in polarization 
from a protoplanetary disk not only depends on the intrinsic de- 
gree of polarization caused by the grains but is also influenced by 
the depolarizing effects of the PSF of the central star. How strong 
this depolarizing effect is depends on the intensity of light scat- 
tered off the disk. Below we discuss the effects of the different 
grain types on the intrinsic polarization image, as well as how 
these different models are affected by telescope effects such as 
diffraction, seeing, and noise. 

In Figs. [T] and [8] we show the different disk and dust types as 
they would be observed with ExPo at the WHT under good see- 
ing conditions, i.e. a seeing at 500 nm of 0.8". There is a clear de- 



pendence of the detectibility on the particle size distribution and 
the vertical disk geometry. The larger and the more settled disks 
are much harder to detect, since they are intrinsically fainter. 

The polarized intensity is the product of the total intensity 
and the degree of polarization. Due to the opening angle of the 
disk, the back side of the disk is closer to 90° scattering angle 
and thus has a higher degree of polarization (see also Fig. [2]). 
This can be seen in the high-resolution images by the longer 
vectors in the upper parts of the images. The front side of the 
disk is, however, much more forward scattering and thus has a 
higher total intensity. Because the central star is the dominant 
intensity source in the observational images, the degree of po- 
larization now is basically the polarized intensity divided by the 
total intensity of the central star and not by the total intensity 
of scattered light. Looking at the observed settled disk models 
(i.e. Fig. [8|, it can be seen that in the images for the compact 
grains, the degree of polarization still dominates (also see the 
length of the polarization vectors), and the back side of the disk 
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Figure?. Same as Fig.|4j but simulated as observed with ExPo at the WHT. Color scale is logarithmic, in arbitrary flux units, and 
the same in all images. 
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Figure 8. Same as Fig.ji] but simulated as observed with ExPo at the WHT. 
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is somewhat brighter in polarized intensity. Due to the highly 
forward- scattering nature of the flufl'y aggregates, in that case 
the difl'erences in total intensity dominate, and the front side of 
the disk is brighter in polarized intensity and in degree of polar- 
ization. 

Since we run the fu l l Ex Po data analysis pipeline as dis- 
cussed in|Canovas et al.|(|2011), we also correct for any instru- 



mental polarization. This is done by assuming that the signal 
from the central resolution element, i.e. the central star, should 
be unpolarized. In the simulation of the instrument we do not add 



polarization caused by the telescope, but we include instrumen- 
tal polarization due to difl'erent transmission coefl&cients (see 
Appendix). Nevertheless, the central resolution element is polar- 
ized because of scattering from the innermost regions of the disk. 
Since the disk is inclined, scattering by this innermost region has 
a non-zero, integrated polarization (see also Fig.fTOl). This polar- 
ization is then 'corrected' under the assumption that it is instru- 
mental, and thus subtracted from the entire image. This causes 
the vertical depolarized region in the images where the outer disk 
is weak. This also causes the well oriented polarization vectors 



M. Min et al.: The effects of disk and dust structure on observed polarimetric images of protoplanetary disks 



9 



in the noise regions of the images for the cases where the disk is 
very weak (especially clear in Fig.[8|. The amplitude of this ef- 
fect strongly depends on the seeing conditions. In Fig. [9] we show 
two disk models under different seeing conditions. It is clear that 
the appearance is different and the overcorrection of the instru- 
mental and inner disk polarization is larger in the case of worse 
seeing. 



In Fig. 10 we show the degree of polarization of the inner 
resolution elements as a function of the size of the resolution 
element. We do this by integrating both Stokes I and Q out to 
a distance R from the central star. The degree of polarization 
shown in Fig.fTOlis the ratio of these two integrated values. This 
shows that increasing the spatial resolution also significantly re- 
duces the effect of polarization from the central resolution ele- 
ment. For images limited by an atmospheric seeing of ~ T', we 
see that the inner resolution element can be polarized by up to 
-0.8%. By adding an extreme adaptive optics system increasing 
the spatial resolution to 0.02'', this degree of polarization can be 
reduced significantly. In the case of fluffy aggregates that have 
grown to significant sizes, we see that the integrated degree of 
polarization is very low. This is because the total scattered in- 
tensity is dominated by forward scattering, which has a very low 
degree of polarization. Although the polarization of the central 
resolution element is low, the spatially resolved degree of polar- 
ization is also extremely low. The large variations in the curves 
in Fig. 10 show that we do have a good diagnostic for the grain 
and disk properties when we have sufficient resolution and the 
images are properly calibrated. 



4.1. Detectability 



For a disk to be detectable with an imaging polarimeter at a 4m 
class telescope without adaptive optics, we need the disk to be 
extended enough to escape the PSF of the central star and bright 
enough to rise above the photon and readout noise. Treating the 
models as observations we have determined which disk and dust 
parameters allow detection with this instrument. For this we as- 
sumed a total integration time of about five minutes per polariza- 
tion state. The visual magnitude of the star, as derived from the 
input parameters in Table[2j is V = 7. All disks are assumed to be 
at a distance of 150 pc. The detectability of the disks is based on 
the signal-to-noise ratio (S/N) of the polarized intensity images. 
The orientation of the polarization vectors depends on this pa- 
rameter, so the polarization vectors calculated for an image with 
poor S/N will be randomly oriented (see, for instance. Fig. [7] 
bottom righthand image). On the other hand, an image with high 
S/N will show a clear pattern on the orientation of its polariza- 
tion vectors. For disks that are hardly detected, we also looked 
at the images for Stokes Q and U to see if a so-called 'butterfly' 
pattern, typical of polarimetric images of circumstellar material, 
emerges above the noise. 

The results are summarized in Tabled We see that disks that 
have undergone significant grain growth are hard to detect, es- 
pecially when the turbulence in the disk is not strong enough 
to mix up the larger grains to the disk surface layers, i.e. in the 
settled disk models. Also, we notice that the seeing conditions 
play a major role in the detectability. By eliminating this fac- 
tor, an adaptive optics system will greatly increase the detection 
capabilities of imaging polarimeters. 



^max 


0.3 jum 


30/um 3000 yL/m 


Good seeing (ro = 10 cm) 


Flared disk, compact grains 


+ 


+ + 


Flared disk, fluffy aggregates 


+ 


+/- 


Settled disk, compact grains 


+ 


+ 


Settled disk, fluffy aggregates 


+ 


+ 


Bad seeing (ro = 5 cm) 



Flared disk, compact grains + + 

Flared disk, fluffy aggregates + 

Settled disk, compact grains + +/- 

Settled disk, fluffy aggregates + 

Table 3. Detectibility of the different disk and dust types with 
ExPo at the WHT. See Table|2]for the details of the model. 



ro = 10 cm (good seeing) ro = 5 cm (bad seeing) 
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Figure 9. Comparison of the polarized intensity maps as ob- 
tained under different seeing conditions. Upper panels are for 
a flared disk, lower panels for a settled disk. All images use 
<^max = 30//m and compact grains. 



4.2. Disk gaps 

Many systems have been observed, directly or indirec tly, to have 
depleted inner disks or cleared out gaps (see e.g. 



Bouwman 

et al. 2003[ [Brown et al.|p07 i Espaillat et al.|2008; Thalmann 



et al 



et al. 



20101 IMerin et al.||2010tlAndrews et al.||2011t IVerhoeff 
201 One of the possible causes for this is clearing by a 



planet in the disk, and thus direct detection of such disk gaps is 
very important for constraining planet formation scenarios and 
timescales. To see if we can detect such disk gaps and see what 
it would generally look like in a polarimetric imaging observa- 
tion, we ran a disk model with a cleared out gap from 120 to 
1 80 AU and an inner disk with a depletion of a factor of ten for 
radii smaller than 120 AU (a parameterization very similar to the 
one adopted by Andrews et al.||2011) . We did not aim to span 
parameter space for this case, but rather to provide an example 
of inner disk clearing. The radii of the gap were chosen to be 
close to what one would expect to be able to detect given the 
spatial resolution we had. Tests with closer in gaps show that, 
indeed, the radii presented here is at the edge of the detection 
capabilities given the simulated spatial resolution. 

The images for compact grains with ^max = 30 yum are shown 
in Fig. [TT] where we show both the full-resolution images (left) 
and the simulated observations (right). For the instrument simu- 
lations we increased the integration time by a factor of ten to sup- 
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Figure 10. The integrated degree of polarization as a function of 
the radius of integration. Upper panel is for the compact grains, 
lower panel for the fluffy aggregates. Black curves are for the 
flared disk models, grey curves for the settled disk models. 



press the photon noise in the inner disk regions. As we can see, 
the depletion of the inner region is clearly detected. Comparing 
to the lefthand panels of Fig. [9] (which shows the same mod- 
els without depleted inner disks), we see a very clear difference. 
We do note that the inner disk is undetectable for the flaring 
disk case in the simulated observations, while it is still clearly 
visible in the full-resolution images. This is caused by the fact 
that the entire inner disk is inside the central resolution element 
(0.8" seeing). As a consequence the polarization signal caused 
by the inner disk is considered instrumental by the data reduc- 
tion pipeline and is subtracted. We also computed images where 
we placed the gap closer in, and we find that an inner radius of 
the outer disk of ~ 180 AU (or 1.2") is close to the limit where 
we are still able to clearly identify the disk gap in the image. For 
radii of the inner edge of the outer disk larger than ~ 100 AU, 
we do see significant changes to the image, but they become in- 
creasingly more difficult to identify as a clear signal from a disk 
gap as the radius decreases. 

5. Discussion and conclusions 

The appearance of circumstellar disks in polarized light depends 
heavily on the properties of both the disk and the composing 
dust particles. We modeled the observable signature of proto- 
planetary disks surrounding young, intermediate-mass stars with 
various dust and disk properties. We focus on the effects of grain 
structure and size-dependent grain settling. We also studied the 
effects of atmospheric blurring and observational data process- 
ing. We conclude the following. 

- There is a significant difference between the appearance of 
a disk with only compact grains and a disk containing fluffy 
aggregated dust particles. 



full resolution 



ro = 1 cnri (good seeing) 



|P=40% 




P=40% 




P= 4% 



P=1 .0% 



Figure 11. Full resolution images (left panels) and simulated ob- 
servations (right panels) of models with a depleted inner disk. 
The innermost regions (radii < 120 AU) are depleted by a fac- 
tor of 10, while from 120 to 180 AU the disk is cleared of dust. 
Upper panels are for a flared disk, lower panels for a settled disk. 
All images are using ^max = 30yum and compact grains. The 
color scale is altered with respect to the previous images to en- 
hance details in the bright inner disk region. 



- The commonly employed homogeneous sphere model gives 
significantly different results from more realistic models, be- 
cause the observed signal depends strongly on the properties 
of these more realistic models, thus providing a potentially 
powerful diagnostic. The proper use of this diagnostic is not 
always trivial. 

- Due to the extreme forward-scattering nature of fluffy aggre- 
gates, where only a small fraction of the radiation is scattered 
towards the observer, a disk containing these kinds of parti- 
cles appears significantly dimmer in scattered light. 

- The diagnostic power of imaging polarimetry is a strong 
function of spatial resolution. Atmospheric seeing blurs 
many details needed for a thorough analysis of dust prop- 
erties, and thus mainly allows the global disk parameters to 
be deduced. 

- The 'ro und about effect' for edge-on disks as discussed by 
"Murakawa ( 2010) for near-IR polarimetry was only found in 
our study for an extremely flaring disk containing only small 
grains. Since the effect appears in the very low-intensity re- 
gions of the image, it disappears when we consider observa- 
tions with limited spatial resolution. 

- The accuracy of the correction for instrumental polarization 
using the central star as calibration source strongly depends 
on the spatial resolution. In the seeing-limited case, the cen- 
tral resolution element can have a degree of polarization of 
up to ~ 0.8%. This degree of polarization quickly decreases 
when the spatial resolution of the observations is increased. 



Appendix A: Instrument simulation 

Here we discuss the details of the instrument simulations. The 
instrument simulator was developed to simulate observations 
with the Extreme Polarimeter (ExPo), which is a visitor instru- 
ment of the 4.2 m WHT. 
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A.1. Wavefront generation 

To generate wave fronts, we make use of a model of turbulence 
first proposed b y |Kolmogorov (|1941 ) and later developed by 
Tatarskil (|1961|) andlFried (1966|). According to this model, a 
flat wavefront traveling through the turbulent atmosphere will 
modify its phase, while the change in its amplitude is negligible 
compared to the phase fluctuations. After it travels through the 
telescope, the amplitude of the wavefront is proportional to the 
telescope's aperture. The power spectrum of the phase fluctua- 
tions of a wavefront traveling through the atmosphere is given 
by 

= ^r"/^ (A.1) 



where ro is Fried' s Parameter and k the wave number. The 
parameter tq represents the circular aperture over which the 
wavefront phase variance is equal to Irad^ for the case of 
Kolmogorov turbulence. Therefore, the higher ro, the more sta- 
ble the atmosphere, or, in other words, the better the seeing is. 
For instance, a value of ro of about 10 cm produces 1 arcsec 
seeing at A = 500 nm (i.e. [Nightingale & Buscher||199i| ). The 
value of ro depend s on the wavelength, following a power law: 



^0 



^6/5 piguj-e 



A.l 



shows an example of an uncalibrated 
wavefront phase generated from Eq. |A.1[ 

The timescale at which one of these wavefronts remains con- 
stant, the coherence time, at optical wavelengths is only a few 
milliseconds ( Vemin & Munoz-Tunon| 1 994} |Law|2006| |Garcia-| 
Lorenzo et al._2009) . The single frame PSF has two main com- 
ponents: 

- Tilt, producing random motion of the whole image. 

- Roughness, producing the observed speckle pattern. 

When the exposure time is similar to the coherence time, the 
image motion due to the tilt component can be removed in the 
data processing by proper centering of each frame. Noll ( 1916) 
showed that the variance cr^ of the wavefronts can be expressed 
in terms of the telescope diameter, D, and the Fried parameter, 

^0 



=0.134 



(A.2) 



The standard deviation a of the wavefront phase if, expressed 
in wavelength units, can be obtained from the variance as ctx = 

0.134i/2(Z)/ro)5/6 '^^^'^ 



cta 



In 



0.0582 







(A.3) 



In our simulation, D is the diameter of the WHT main mir- 
ror (4.2 meters), and three difl'erent values of ro, (5, 7, and 10 
cm) are used to account for bad, medium, and good seeing at a 
reference wavelength of 500 nm, /I500. From Eq |A.2| and the de- 
pendence of ro with wavelength, (Jx can be computed at difl'erent 
wavelengths from our reference as 



^500 



(A.4) 



W avefro nts at di fferent wavelengths are calibrated according 
to Eqs. A.2 and A.4 A set of wavefronts at difl'erent wavelengths 
is then generated for each of the three difl'erent ro values used in 
this simulation. 



Parameter 


Value 


Telescope diameter 


4.2 [m] 


Wavelength range 


450-700 [nm] 


Exposure time 


0.028 [seconds] 


Pixel size 


0.078 ["/pixel] 


Grid size 


2048 X 2048 [pixels] 



Table A.l. Simulation parameters 



A.2. PSF generation 

To generate a PSF, we first simulate the aperture of the WHT, 
taking the spiders of the telescope and the central obscuration of 
the main mirror into account. Once the telescope aperture and 
the wavefront are computed, a monochromatic PSF is calculated 
according to Eq. [2 A broadband PSF (PSF/^/,) is calculated as 
the sum of monochromatic PSFs calculated over the range of 
400 nm to 700 nm, in steps of 10 nm: 



700 



PSF, 



(A.5) 



i=400,10 



This PSF produces the speckle pattern obtained when observ- 
ing at exposure times that are close to the coherence time (a few 
milliseconds). The PSF measured when observing at longer ex- 
posure times will be the sum of the short-exposure broadband 
PSFs. To simulate ExPo-like PSFs, we set the coherence time to 
9.3 milliseconds, and then generate a 28 millisecond exposure 
PSF as the sum of three statistically independent short-exposure 
broadband PSFs: 



PSF 



ExPo 



(A.6) 



These calculations are computed on a 2048 x 2048 pixels grid. 
The pixel size of this simulation is determined by the Nyquist 
frequency: Ny = (A/ 2D) • 206265, which produces a pixel size 
of 0.0122 ["/pixel], for a 4.2 meter telescope and a central wave- 
length of 500 nm. The original 2048 x 2048 pixel grid is then 
binned to an smaller grid to produce si mulate d images at the 
ExPo pixel size (0.078 ["/pixel]). Figure A.2 shows an exam- 
ple of a PSF for both simulated (left) and real (right) data (0.8" 
seeing, 0.028 seconds exposure time). 

A total of 100 diff'erent PSFexPo are produced for each of the 
three diff'erent seeing conditions tested here. Figure A.3 shows 
an example of three diff'erent broadband PSF's calculated for bad 
(left), normal (center), and good (right) seeing. 

A.3. Dual beam simulation 

As a dual beam instrument, ExPo produces two simultaneous 
images with opposite polarization states that are imaged onto a 
CCD. We refer to these two images as left and right images. 
These measurements are modulated by a Ferro-electric liquid 
crystal (FLC), which switches between two orthogonal polar- 
ization states (A and B). Two diff'erent frames, containing four 
diff'erent images are produced at the end of one FLC cycle: 



Al = 0.5(7 + P)o*PSFa, 

Ar = 0.5(7 -P)o*PSFa, 
Bl = 0.5(7 -P)o*PSFb, 
Br = 0.5(7 -hP)o>^PSFB, 



(A.7) 
(A.8) 
(A.9) 
(A. 10) 
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Figure A.l. Phase of a wavefront following the Kolmogorv 
model of turbulence. 
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Figure A.2. Left: Simulated PSF with seeing ^ 0.8", 0.028 sec- 
onds exposure time. Right: Observed PSF under the same con- 
ditions. 
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Figure A.3. Simulated speckle pattern for different values of tq, 
given in centimeters. 



where neither instrumental polarization nor instrumental effects 
are considered for the sake of simplicity, and /q and represent 
the total intensity and polarization (Stokes Q or Stokes U) of 
the observed target, as seeing without atmosphere and diffraction 
effects, respectively. PSFa and PSFb are the short exposure PSF 
for the A and B frames, respectively. These images are given 
in mJy units, and they are converted to counts by considering 
the telescope area, exposure time, filter transmission (Johnson 



Coefficient Value 



Tl 


0.48 


Tr 


0.52 


Ta 


0.998 


Tb 


1.002 



Table A.2. Transmission coefficients measured from ExPo in- 
strument. 



V filter is used here), atmospheric + instrument absorption, and 
CCD efl&ciency. 

The PoiDEV function from the NASA IDL library is used to 
simulate the photon noise for each of these images. Real readout 
noise from ExPo measurements is added to the simulated data. 
Instrumental polarization is simulated by including two trans- 
mission coefl&cients for each of the beams (Tl,Tr) and another 
two coefficients to account for the FLC transmission (Ta,Tb). 
The value of these coeflftcients is listed in Table IA.2I To mimic 
guiding problems, a random shift with a maximum amplitude of 
ten pixels is applied to each of the simulated images. The final 
simulated images are then described by 

Al = Ta • Tl • (/ + P)o * PSFa • M(x, y) + PHal + ROal, (A. 1 1) 
Ar = Ta • Tr . (/ - P)o * PSFa • M(x, y) + PHar + ROar, (A. 12) 
Bl = Tb • Tl • (/ - P)o * PSFb • M(x, y) + PHbl + RObl. (A. 13) 
Br = Tb • Tr • (/ -h P)o ^ PSFb • M(x, y) -h PhsR + RObr, (A. 14) 

where M(x, y) represents the image- shifting function, Ph is the 
photon noise, and RO represents the readout noise. 

We finally run the full data reduction pipeline for the simu- 
lated observations to obtain the final images presented in Figs. [7] 
[8land|9] 
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